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ABSTRACT 

By exploiting the exceptional high-resolution capabilities of the near-IR cam¬ 
era GSAOI combined with the multi-conjugate adaptive optics system GeMS at 
the GEMINI South Telescope, we investigated the structural and physical prop¬ 
erties of the heavily obscured globular cluster Lillcr 1 in the Galactic bulge. We 
have obtained the deepest and most accurate color-magnitude diagram published 
so far for this cluster, reaching K s ~ 19 (below the main sequence turn-off level). 
We used these data to re-determinc the center of gravity of the system, finding 
that it is located about 2.2" south-east from the literature value. We also built 
new star density and surface brightness profiles for the cluster, and re-derived 
its main structural and physical parameters (scale radii, concentration parame¬ 
ter, central mass density, total mass). We find that Liller 1 is significantly less 
concentrated (concentration parameter c = 1.74) and less extended (tidal ra¬ 
dius r t = 298" and core radius r c = 5.39") than previously thought. By using 
these newly determined structural parameters we estimated the mass of Liller 1 
M tot = 2.3 ^3 x 1O 6 M 0 (M tot = 1.5 x 1O 6 M 0 for a Kroupa IMF), which is 
comparable to that of the most massive clusters in the Galaxy (oj Centauri and 
Terzan 5). Also Liller 1 has the second highest collision rate (after Terzan 5) 
among all star clusters in the Galaxy, thus confirming that it is an ideal environ¬ 
ment for the formation of collisional objects (such as millisecond pulsars). 

Subject headings: Globular Clusters: Individual (Liller 1) - Stars: evolution, 
Technique: photometry 
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1. Introduction 


Liller 1 is one of the 7 globular clusters (GCs) located within 1 kpc from the Galactic 
centre, at a Galactocentric distance of only 0.8 kpc, and very close to the galactic plane 
(Z = 354.84, b = —0.16; Harris 1996 


, 2010 edition, hereafter H96). This region of the 


Galaxy is strongly affected by a large foreground extinction: the average color excess 


E(B — V) estimated by c 


3.00 to 3.09 (jFrogel et al 


ifferent authors in the c 


(119951! : 


Qrtolani et ah 


significant evidence of differential reddening (Ortolani et al. 


irect i on of t he cluster ranges fro m about 


1996 


2001 


Valenti et ah 


(12 0101) with 


2001). Such a large value of 


the extinction has severely hampered the observations of Liller 1 in the optical bands. 
Indeed very few papers about this cluster can be found in the literature: the only_available 


optical color-magnitude diagram (CMD) of Liller 1 is the one published by 


Ortolani et ah 


( 19961) by using ESO-NTT data. V 
is the (K,J-K) CMD published by 


ore m genera 


Frogel et al. 


^the first CMD ever obtained for Liller 1 


(119951 !. However in both cases data are 


very shallow and the authors were able to sample only the brightest part of the red giant 


branch 


(RGB). The 


Davidgel (2000) and 


deepest photometry of Liller 1 published so far, was obtained by 


Ortolani et al. 


(1200111 using the Canada-France-Hawaii Telescope and 


the Hubble Space Telescope, respectively. However, these samples are not deep enough to 
properly characterize the main sequence turn-off (MS-TO) region of the cluster. A more 


recent photometric analysis of Liller 1 was performed by 


Valenti et al. 


(120101! by using 


the large near-IR wide held of view (FoV) imager Sofl mounted at the ESO-NTT. The 


1 Based on observations obtained at the Gemini Observatory, which is operated by the Association of 
Universities for Research in Astronomy, Inc., under a cooperative agreement with the NSF on behalf of 
the Gemini partnership: the National Science Foundation (United States), the National Research Council 
(Canada), CONICYT (Chile), the Australian Research Council (Australia), Ministerio da Ciencia, Tecnologia 
e Inovagao (Brazil) and Ministerio de Ciencia, Tecnologia e Innovation Productiva (Argentina). Based on 
observations gathered with ESO-VISTA telescope (program ID 179.B-2002). 
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obtained CMD samples the brightest portion of the RGB reaching the cluster red clump 
level (K s ~ 14). 

Liller 1 is one of the most metal-rich GC in th e Galac tic bulge. In fact b y using 


high-resolution IR 


spectroscopy ([Origlia et al. 


1 99 7). both 


Origlia et al. 


(120021 1 and 


Stephens fe Fro sell (1200411 measured a metallicity of about one h alf solar ([Fe/Hl = — 0.3 


dex). A similar value ([Fe/H]= —0.36 dex) was also estimated by 


Ferraro et al. 


Valenti et ah 


(2010) from 


2000) characterizing the location and the 


a set of photometric indices (see 
morphology of the RGB in the (K x , J — K s ) CMD, and calibrated as a function of the 


metallicity (Valenti et al. 


2004a,b). 


One of the most astonishing ch aracte ristics of Liller 1 is the extremely large value of 


the collision rate parameter. 


Verbunt fe Hut 


(119871) showed t hat Liller 1 has the s econd 


highest value of stellar encounter rate (after Terzan 5; see also 


Lanzoni et al. 


2010) among 


all star clusters in the Galaxy, thus suggesting that it represents an ideal environment where 
exotic objects, generated by collisions, can form. In fact, it is commonly believed that 
dynamical interactions in GCs facilitate the formation of close binary systems and exotic 
objects like Cataclysmic Variables (CVs), Low Mass X- 

er Stars (BSSs) (e.g. 


Pulsars (MS’ 

V) an< 

Blue 

Stragg 

Ferraro et al 

2001, 

2009a, 

2012; 


(see also 


Ransom et al 


2005 


C-rav Binaries 

XMXBs), Mi 

isecond 

Bailvn 

1992 


3 aresce et al. 

1992 


Poolcyi & Hut 

2006 

). Moreover, 



Hui et al .l (120101) found t hat clusters with large collisional parameters and high metallicity 


Bcllazzini et al. 


19951) usually host more MSPs. Indeed T erzan 50 hosts the largest 


population of MSPs among all Galactic GCs (Ransom et al 


2005). A strong y-ray emission 


has been recently detected in the direction of Liller 1 by the Large Area Telescope (LAT) 


2 Note that this stellar system is suspected to not be a genuine GC, because it liar- 


bors at least three stellar populations with d ifferent iron abundances (IFerraro et al. 


Orielia et al. 

2011, 

2013; 

Massari et al. 

2014) 


2009b 
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on board the Fermi Telescope (Tam et al. 


2 01ll) . This is the most intense emission detected 


so far from a Galactic GC, again suggesting the presence of a large number of MSPs. 
Ho weve r, no dir ect r adio detection of these objects has been obtained so far in this system 


(IRansom et ah 


2005 b The only exotic object identified in the cluster is the rapid burster 


MXB 1730-335, an LMXB observed to emit radio waves and type 1 and type II X-ray 


bursts (IHoffman et al. 


19781 ). It seems to be located in the c entral region of Li 


optical/IR counterpart of this object has been found so far fjHomer et al. 


2001 ). 


ler 1, but no 


In this paper we present the deepest near-IR photometry of Lillcr 1 yet obtained and 
new determinations of its main structural and physical parameters. In Section [21 we discuss 
the observations and data analysis. Section [3] presents the results obtained, including new 
determinations of: the distance modulus and reddening (Sect. 13.11) . the center of gravity 
(Sect. 13.21) . the star density profile, the structural parameters and the surface brightness 
(SB) profile (Sect. 13. 3p . the total luminosity, total mass and central mass density of the 
cluster (Sect. 13.4[) . The summary and conclusions are presented in Section [4] 


2. Observations and data analysis 


The main photometric data used in the present study consist of a set of high-resolution 
images obtained with the IR camera Gemini South Adaptive Optics Imager (GSAOI) 
assisted by the Gemini Multi-Conjugate Adaptive Optics System (GeMS) mounted at the 
8 m Gemini South Telescope (Chile), in April 2013 (Program ID: GS-2013-Q-23; PI: D. 
Geisler). GSAOI is equipped with a 2 x 2 mosaic of Rockwell HAWAII-2RG 2048 x 2048 


pixels arrays. 


flNeichel et al. 


t coy ers a global FoV of 85" x 85" on the sky, with a resolution of 0.02"/pixel 


2014). We sampled the central region of Liller 1 with a mosaic of multiple 


exposures acquired with a dithering pattern of a few arcseconds. Seven and ten exposures 
have been acquired in the J and K s bands, respectively, with an exposure time t exp = 30 s 
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each. The entire data set was acquired in an interval of time of ~ 90 mins. A constellation 
of five laser guide stars combined with three natural guide stars has been adopted to 
compensate the distortions due to the turbulence of Earth’s atmosphere. We have verified 
that the Full Width at Half Maximum (FWHM) is quite stable allover the FoV with a 


maximum variation of about 10 — 15%, in agreement with what found by 


Neichel et al. 


(2014). The average FWHM (i.e. estimated on the entire FoV) varies from 0.09" to 0.13'' 


for the images acquired in the J band and from 0.08" to 0.12" for the K s ones. These values 
are only slightly larger than the diffraction limit of the telescope. In Figure [Tj we show a 
two-color image of Lillcr 1, obtained by combining GEMINI J and K s band observations. 

As discussed in Section [3] below, in order to derive the radial star density and SB 
profiles of Liller 1, we complemented this data set with a near-IR catalog from the VISTA 
Variables in the Via Lactea Survey (hereafter VVV, 


Minniti ct al. 


2010 ; 


2011) complemented with the Two Micron All Sky Survey (2MASS; 


Catelan et ah 


Outri et al 


2003) 




For the GEMINI data set, a standard pre-reduction procedure, using IRARj tools, 
was applied to the raw images to correct for flat fields and bias, and to perform the 
sky-subtraction. The photometric reduction was then carried out via point-spread function 
(PSF) fitting technique in each chip of each image independently, using DAOPHOT 


((Stetson 


19871) . a package developed to perform accurate photometry in crowded fields. 


The PSF has been modeled by selecting about one hundred bright and isolated stars 
uniformly distributed in each chip, and by using the DAOPHOT/PSF routine. We allowed 
the PSF to vary within each chip following a cubic polynomial spatial variation. The 
best-fit PSF analytic models obtained for the J and K s images, selected on the basis of a 


3 IRAF is distributed by the National Optical Astronomy Observatory, which is oper¬ 
ated by the Association of Universities for Research in Astronomy, Inc., under cooperative 
agreement with the National Science Foundation. 
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y 2 te st, are a Penny function ( Penny 1976) and a Moffat function with (3 = 1.5 ( Moff at 
19691 ). respectively. The PSF models thus obtained were then applied to all the star-like 
sources detected at a 3a level from the local background by using ALLSTAR. We have 
derived in this way the stellar instrumental magnitudes. Then, starting from the star lists 
thus obtained and to fill the gaps among the GSAOI chips, we created a master star list 
containing all the stars measured in at least one J and K s image. 


As done in other works (see, e.g., 


Dalessandro et al. 


2014. and r eference there in), the 


master-lists thus created have been used as input for ALLFRAME (IStetso n 


19940- The 


files obtained as output have been combined to have a complete catalog with the J and 
K s magnitudes and the positions of all detected stars. For every stellar source, different 
magnitude estimates have been homogenized and their mean values and standard deviations 


have been finally adopted as 

(see 

Ferraro et al. 

1991, 

1992) 


The distribution of the photometric errors as a function of the J and K s magnitudes 
are shown in Figured Errors vary from 0.01 - 0.02 for bright RGB stars, and they reach 
values a rm K 0.1 cind. (J rrn j 0.3. As expected, given the overall stability of the PSF, we 
do not observe any significant radial trend. 


The instrumental magnitudes have been converted into the 2MASS photometric system 
by using the stars in common with the publicly available catalog obtained with SOFT 


(Valenti et ah 


2010 ). To minimize the effect of blending, photometric errors and saturation, 


we used only stars at a distance larger than 30" from the cluster center and with 2.8 < 
(SOFI K s - GeMS K s ) < 3.6 and 2.6 < (SOFI J - GeMS J) < 3.6, respectively. The J 
and K s calibration curves are shown in Figure [3j As can be seen, we used an iterative 
sigma-clipping procedure to estimate the median of the (SOFI K s - GeMS K s ) and (SOFI 
J - GeMS J ) distributions for the stars satisfying the criteria listed above. The r.m.s of the 




















best-fit relation is 0.1 and 0.08 in K s and J, respectively. The values thus obtained (Figure 
ED have been applied to the instrumental GeMS magnitudes. The same stars have been 
used also to put the instrumental positions onto the absolute astrometric system. 


We retrieved VVV images obtained in the direction of Lillcr 1 from the Vista Science 


Archive website 4 5 


The data set is composed of one exposure per filter and it covers a 


total FoV of 3 3 A x 33A. Data were prereduced at the Cambridge Astronomical Survey 


Unit (CASU 


fl 


with the VIRCAM pipeline (Irwin et al. 


20041 1 . We then p erformed 


PSF-fitting photometry by using the DAOPHOT based VVV-SkZ_pipcline (jMauro et al. 


2 01311 o n the single 2048 x 2048 pixels chip extracted from the stacked VVV pawprints 


flSaito et al. 


2012). A quadratically variable Moffat function with /3 = 3.5 has been adopted. 


The VVV magnitudes and instrumental coordinates were then 
photometric and astrom etric system (see details in 


2011 


Mauro et al. 


2013). 


s were then rep or 

ed to the 2MASS 

Chcne et al. 

2012: 

Moni Bidin et al. 


We estimated the differential reddening in the direction of Liller 1 by adopting a 


procedure similar to that discussed in 


Massari et al. 


(2012), with the only difference that 


here we used the RGB instead of the MS stars as reference sequence. In particular, we 
divided the GSAOI FoV in a grid of 6 x 6 cells, each 18"-wide. We considered RGB 
stars approximately in the magnitude range 14.5 < K s < 16.5: the upper and the lower 
thresholds of the selection box have been set running parallel to the reddening vector. 
Considering the RGB stars within each cell, we estimated the median color. The cell with 
the bluest color (nominally corresponding to the lowest extinction value) has been adopted 
as reference. The relative color shift h[( J — K s )\i of each zth cell, is then defined as the 
shift needed to make the median color of the ith cell match the reference cell color. Then, 


4 http://horus.roe.ac.uk/vsa/ 

5 http://casu.ast.cam.ac.uk/ 
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from the value of S[(J — K s )]i we derived the corresp onding 5\E(B — V)] 


the extinction coefficients Rj = 0.87 and Rk s = 0.35 flCardelli et al.l 1989). Extinction 


by adopting 


variations as large as 5[E(B — V)\ = 0.34 mag have been measured in the direction of the 
cluster. 


Figure H] shows the differential reddening corrected ( K S1 J — K s ) CMD of Liller 1. 
This is the deepest and most accurate CMD ever obtained for this stellar system, reaching 
K s ~ 19 and thus sampling the MS-TO region. Its well defined red clump, clearly visible 
at K s ~ 14.2, is a typical feature of a metal rich GC. Unfortunately, stars lying along 
the brightest portion of the RGB (K s < 12) are saturated in all the available images. As 
apparent from Figure [5l Liller 1 suffers from significant held contamination mainly from the 
Galactic bulge and disk stars, which define a blue plume clearly visible in the bluest portion 
of the CMD for (J — K s ) < 1.5. While the features of the CMD will be discussed in detail 
in a forthcoming paper (hereafter, paper II), here we focus on the derivation of the star 
density profile of Liller 1 and the estimate of its main structural and physical parameters. 


Results 


3.1. Distance and reddening of Liller 1 


The most recent determinations of the distance and reddening of Liller 1 (Valenti et al. 


2010) have been obtained on the basis of a relatively shallow ( K s , J — Jl s )-CMD, reaching 
only the red clump level. Since the GEMINI photometry presented here is sig nificantl y 


neper, we used 


(120071. see also 


our new data to re-deter mine both these quantities. Following 


Valenti et ah 


Dalessandro et ah 


2008a ). we used a differential method, consisting in 


the comparison between the CMD and the luminosity function of Li 


ler 1 with those o 


2002 


a reference cluster, NGC 6553, a bulge GC with similar metallicity (jOriglia et ah 
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Melendez et al. 

200 

3; 

Alves-Brito et al. 

2006 

NGC 6553 is from 

Valenti et al. 

(2007, 

see a 


Valenti et al 


Ferraro et al. 


2007). The data set used for 


20001 ). We found that, 


in 


order to align the RGB mean ridge line and the magnitude level of the red clump of the 
two clusters, a color shift 5(J — K s ) = 1.2 and a magnitude shift 5K S = 1.9 must be applied 
to the sequences of NGC 6553 (see Figure [6]). Thus, by adopting [m — M) 0 = 13.46 and 


E(B — V) = 0.84 for NGC 6553 ( Valenti et ah 2010), we obtain (m — M)k = 15.65 ± 0.15 
and E(B — V) = 3.14 ± 0.20 for Lillcr 1. The relative uncertanties on (m — M)k and 
E(B — V) are obtained by taking into account the photometric errors at the magnitude 
level of the red clumps and RGB bumps of Lillcr 1 and NGC 6553, the errors due to the 
photometric calibration (Section [2]) and the histogram bin size. 


It is important to note that the adopted reddening value has been obtained by using 
the differential reddening-corrected CMD of Liller 1, hence it corresponds to the least 
extincted region in the GEMINI FoV (see Section [2]). By averaging the values of E(B-V) 
of each cell considered in Section [2] weighted by the number of reference stars sampled, we 
obtain a mean extinction of E(B — V) = 3.30 ± 0.20 in the direction of Lillcr 1. This yields 
a true (unreddened) distance modulus of (m — M) 0 = 14.55 ± 0.25, which corresponds to a 
distance of 8.1 ± 1.0 kpc, in agreement with previous determinations (Harris 1996; Valenti 
et al. 2010; Ortolani et al. 2007). 


Interestingly. 


the RGB bump (Fusi Pecci et al. 


Figure [6] also s 


rows t 


1990; 


rat the relative posi tions of the red clump and of 


Ferraro et al 


1999a) are quite similar in the two 


clusters. Hence, since the location in magnitude of the RGB bump is quite sensitive to 
metal abundance, such a nice correspondence further supports the evidence that Lillcr 1 
and NGC 6553 share the same chemistry. 
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3.2. Center of gravity 


Using the absolute positions of individual stars in the GEMINI sample, we determined 


the center of grav i ty (C A-av) of L iller 1, by following the iterative procec 


Montegriffo et ah 


(11995 


see also 


Ferraro et al. 


2003 


Lanzoni et ah 


ure described in 


20071) . averaging the 


right ascension (a) 

and declination (<5) va 

ues o: 

Ferraro et al. 

(2003 

, see also 

Lugger et al. 

1995; 


properly selected s tars. As discussed in 


Miocchi et al. 


20131 ). the use of individual 


stars to measure the center of gravity and the star density profile provides the most robust 
approach to derive cluster structural parameters, since the counting of individual stars is 
not affected by the presence of a few bright objects (which can, instead, generate spurious 
luminosity clumping in the measurements based on the observed SB). However star counts 
can be affected by incompleteness. We estimated the effect of completeness by means 


described in 


Dalessandro et al. 

(2011 

, see also 

Bcllazzini et al. 

2012) 


2012). We obtained that 


only stars with 12.5 < K s < 15.5 have a completeness C ~ 100% at any distance from the 
cluster center. We used only these stars to study both the center of gravity and the density 
profile (Section 13.31) . 


We estimated the gravity center by considering different sub-samples of stars: we used 
stars lying within different distances (10", 20" and 30") from a first guess center, and with 
three different cuts in magnitudes (K s =15.5, 15.0 and 14.5). The final value of C grav is 
the average of the different estimates. It turns out to be located at aj 2 ooo = 17 h 33 m 24.56 s , 
<5j2ooo = —33°23'22.4", with an uncertainty of ±0.3" and ±0.2" in a and 5, respectively. 
Our center is ~ 2.2" south-east from that reported by H96 and estimated by using the SB 


distribution. 
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3.3. Star density profile and structural parameters 

By adopting the derived value of C grav and the data sets described in Section [2] we 


to 900"). Fo l lowing the 


procedure described in 


20071 ; IDalessandro et al. 


Ferraro et al. 

(1999b 

, see also 

Lanzoni et al. 


2008bl) . we divided the GEMINI and the VVV samples in several 


concentric annuli centered on C grav , each one split into a variable number of sub-sectors. 
The number of stars lying within each sub-sector was counted, and the star density was 
obtained by dividing these values by the corresponding sub-sector areas. The stellar density 
in each annulus was then computed as the average of the sub-sector densities and the 
standard deviation was adopted as the uncertainty. In order to avoid incompleteness and 
saturation biases, different limiting magnitudes have been adopted for the two data sets: 
12.5 < K s < 15.5 for the GEMINI sample and 12 < K s < 15.5 for the VVV data set. 

By guaranteeing an adequate radial overlap (Ar = 30 — 50") among the two samples, the 
annuli in common between the adjacent data sets have been used to join the two portions 
of the profile. At the end of the procedure, the inner portion (r < 50") is sampled by the 
high-resolution GEMINI data, while the outermost portion is obtained from the VVV data 
set. 

The projected density profile thus obtained is shown in Figure [7] (empty circles). The 
distance from the center associated with each point corresponds to the mid-value of each 
adopted radial bin. The projected density at r > 100" is approximately constant, consistent 
with being due to the Galactic field population alone. To estimate the background and 
foreground Galaxy contamination we therefore used the lcr-clipped average of the five 
outermost points (dotted line in Figure [7]). The field-decontaminated star density profile of 
Liller 1 is marked by the solid circles in the figure. The observ ed d ensity distributions of GCs 


are traditionally described by means of King models (IKing 


3 ), 


1966) , even if deviations from 
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this 


dnd of profile have been found in some cases (see, e.g.. iMcLanghlin &; van der Mare! 


2005 


Moni Bidin et al. 

2011; 

Miocchi et al. 

2013 


the observed star density profile of Lil 
an isotropic, single-mass King m odel flKing 


er 1 and c 


and references therein). To reproduce 
erive its structural parameters we used 


described in 


Miocchi et al. 


19661) . and we followed the procedure fully 
( 20131 ). According to a x 2 test, the density profile of Lillcr 1 can 


be excellently reproduced by a King model with core radius r c = 5.39" + q^, concentration 
parameter c = 1.74^o;}g, half-mass radius ry = 30.and tidal (or limiting) radius 
r t = 298 ";“ 

The structural parameters we obtained are significantly different from those quoted by 
H96 based on the SB profiles in I and R bands (r c = 3.6", c = 2.3 and thus r t = 720"). In 
particular, the newly derived parameters show that Liller 1 is significantly less concentrated 
and less extended than previously thought. A reasonably good agreement_is found 


with the r c values estimated from the an alysis of JHK S images by 
V' _x_ n cV'Jfi 


ry = 5.2" ± 0.5' 1 


and by 


Malkan et al. 


Cohen et ah 


(2007, 


(198(J, r c = 7" ± 2"). A plausible explanation 


for the disagreement between the structural parameters derived from optical and near-IR 
observations could be a differential reddening effect, that may produce a spurious distortion 
of the optical SB profile, as well as low S/N in the optical due to the heavy extinction. 

We derived also the SB profiles of Lillcr 1 directly from the images. To this aim, we 
used the wide-held images of the 2MASS data-set, since the VVV images heavily suffer 
fr om saturation. We produ ced J an d K s SB profiles b y following the procedure described 


m 


Dalessandro et al. 


(2012 


, see also 


Cohen et ah 


2007) and by adopting the same center, 


Tt is worth noticing, however, that the approach used in 


Cohen et ah 


( 2007 ) to fit the 


SB profile is quite different from ours, since the authors derived the central SB and the core 
radius of the cluster after adopting r t and C grav from the online database of H96. Hence, the 
results obtained are not directly comparable. 
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and a setup of annuli and sub-sectors similar to the one used to obtain the density profile. 
Basically, for each annulus we summed the photon counts sampled in each sub-sector and 
we assumed as final SB the average value of the different sub-sectors. In order to reduce 
statistical fluctuations due to the presence of a few bright giants randomly distributed in 
the FoV of the J and K 2MASS images, a threshold of 8000 photon counts was adopted 
for each pixel. The counts have been converted into instrumental magnitudes and then 
calibrated by using appropriate zero-points. The contribution of the Galactic background 
was finally estimated for r > 130", obtaining /x^ ck = 12.3 mag/arcsec 2 and /ij ack = 14.5 
mag/arcsec 2 . The SB profiles thus obtained are shown in Figure |U These profiles are well 
fit by the same King model derived from the projected density distribution, thus confirming 
the accuracy of the structural parameters derived in the present work. The values of the 
central SB measured in the J and K s bands, /ij(0) and /i#(0), are listed in Tabled! together 
with all the relevant parameters derived for the cluster. 


3.4. Cluster physical parameters 


Given the value of the effective radius estimated from the best-fit King model discussed 
in the previous Section (r e = 22.4"75g), we derived the integrated K s luminosity of 
the cluster within r e . By using aperture photometry on the 2MASS images and after 
subtracting the background contribution determined as described in Section 13.31 we 
obtained K s (< r e ) = 4.30+q;]J- Hence, by definition of r e (i.e., the projected radius 
including half the total integrated light), the total integrated magnitude of the cluster 
is K„ = 3.55^q} 9. This value is significantly brighter than previous estimates. In fact 


Cohen et ah 


(120071 ) found K s (< 50") = 4.5, while our estimate over the same cluster region 


is 0.6 mag brighter, K s (< 50") = 3.9. We emphasize, however, that our estima 
measure of the integrated K s magnitude over the 2MASS image, while that of 


te is a direct 


Cohen et al. 






15 


(2007) depends on the assumed cluster parameters (which are different from those derived 
here, as discussed in Section! 


By adopting a color excess E(B — V) = 3.30 ± 0.20, a distance modulus 
(m — M )o = 14.55 ± 0.25 (see Section l3Tj) . and a bolometric correction BCj^ = 2.3 


appropriate for a population of intrinsic color (J — K s ) 0 = 0.75 ( Monte gr iffo et ah 19931 1. we 


estimated a bolometric luminosity of about Aboi(< f e ) = 3.45 + q;45 x 10 5 A o , corresponding 
to a total luminosity Aboi = 6.9+o'g x 1O 5 L 0 . From this value, the total mass of Lillcr 1 


can be estimatec 


2000 


Maraston 


bsL 


assuming a mass-to-light ratio M/L^ 0 \ = 3.35 (iMaraston & Thomas 


2003) appropriate for a Salpeter Initial Mass Function (IMF), thus 


obtaining M tot = 2.3 + q;3 x 10 6 M 0 . We also estimated the mass of Liller 1 assuming 
M/Lboi = 2.19 as obtained by iMaraston (j2003l) for a Kroupa IMF. We obtain in this case 


M tot = l-5-°i x 1O 6 M 0 . 


Note that such a large value o 


( Merritt et al 


1997) and Terzan 5 (Ferraro et ah 


the mass is similar to that extimated for u j Cen 


2009a; 


Lanzoni et al. 


20101 ). and rankes 


Liller 1 in the high-mass tail of the distribution observed for GC-like stellar systems in the 
Milky Way. 

Same values of the total mass are obtained from the total integrated A'-band 


magnitude (see above' 


(IMaraston & Thomas 


and 


2000 


adopting mass-t o-light ratios M/L K = 1.56 and M/A^ = 1-03 


Maraston 


2003) for a Salpeter and a Kroupa IMF, respectively. 


We also used the available data to get a new estimate of the central mass density 
of the cluster . To do so, we followed the same procedure described for K-band data 


m 


Diorgovski (119931 ). However, in order to minimize the effects du e to t he strong and 


differential reddening affecting the system, we re-derived eq. (5) of iDiorgovski (119931 1 in the 


appropriate form for A7,-band photometry. By adopting 3.28 as the A) s -band magnitude of 
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the Sun, we obtain: 


where / n is 


of 


log/ 0 , Ks = 0.4(24.852 -p A -(0)], (1) 

he central projected luminosity density in units of L 0 pc -2 . From eq. (4) 


Diorgovski (1993]) and the structural parameters quoted in Section [3731 we then derived 


the central luminosity density of the system in the K- band. By assuming M/L A = 1.56 
appropriate for a Salpeter IMF, we finally converted it into the central mass density, 
p 0 ~ 7.2^;° x 1O 6 M 0 pc -3 . Using instead a Kroupa IMF which gives M/L K = 1.03, we 
obtain p 0 ^ 4-8+L2 x 1O 6 M 0 pc 3 . These values are between 5 and 7 times larger than 


those obtained from the U-band centra] 
M/L v = 5.3 (from 


luminosity density qu oted in H96, by assuming 


Maraston & Thomas 

2000; 

Maraston 

2003) 


4. Summary and conclusions 

By using GeMS+GSAOI at the Gemini South Telescope, we have obtained the 
deepest CMD of the Galactic GC Liller 1 published so far, properly sampling even the 
innermost regions of the system (except for some saturated stars). This allowed us to 
obtain new estimates of the cluster distance and reddening, which essentially confirm 
previous determinations. Significant differences with respect to the literature values, 
instead, have been found for the cluster center and structural parameters: with respect to 
the values quoted in the H96 catalog, Liller 1 turns out to be located ~ 2.2" south-east, 
to be significantly less concentrated (c = 1.74, instead of c = 2.3), to have a larger core 
radius (r c = 5.39", instead of r c = 3.6"), and thus to be less extended overall ( r t = 298", 
instead of r t = 720"). Also its total mass M tot = 2.3 x 1O 6 M 0 and central mass density 
p 0 = 7.2 x 1O 6 M 0 pc -3 ( M tot = 1.5 x 10 6 M 0 and p 0 = 4.8 x 1O 6 M 0 pc -3 for a Kroupa IMF) 
appear to be a factor of a few larger that those derivable from the U-band parameters listed 
in the H96 catalog. 
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From these new estimates, we also re-determined the collisional parameter of the 
system, which can be expressed as T oc /?q 5 x for virialized and King model systems. We 
find that T is ~ 20 - 40 times larger than the one obtained by using the H96 parame ters. It 


corresponds to about half the collisional parameter of Terzan 5 estimated in 


Larizoni et ah 


(2 0 10), and is much larger than the values derived for other massive GCs for which 


the structural parameters have been re-d etermined consistently wi t 


l what is done here 


(NGC 6388, NGC 6266, 47 Tucanae; see 


Mapelli et ah 


2006 


Dalessandro et al. 


2008a; 


Beccari et al 


2006 


, respectively). 


Hence, our analysis confirms the previous suggestion (Verbunt & Hut 1987) that the 


collisional parameter of Liller 1 is the second highest (after that of Terzan 5) among all 
Galactic GCs. In stellar systems with high values of T, large populations of “collisional” 
objects (like LMXBs, MSPs, BSS and CVs), whose formation is promoted by single/binary 
encounters, should not be unusual (also depending on the availability of the progenitors). 


In this respect, it is interesting to note th at, w 


lile 


of MSPs ever observed in a GC ( Ransom et ah 


erzan 5 hosts the largest population 


20051), no MSPs have been c 


in Liller 1, the only exotic object being in fact a LMXB flHomer et al. 


2001 


etected so far 
, see also the 


Introduction). The lack of detection of large populations of MSPs, LMXBs and CVs might 
be due to observational biases. In fact Liller 1 is located in a region of the Galactic plane 
where diffuse emission is very strong. In this respect, deeper MSP searches, perhaps with 
the Square Kilometer Array, could help understanding the real MSP content of the system. 


On the other hand, the strong y-ray emission (the mos 


detected in the direction of the cluster ( Tam et a_L 


intense among all Galactic GCs) 


2011) may suggest the possible presence 


of numerous, but still hidden, MSPs. Another possibility, emerging from the results 
pr esen ted in this work, is that the y-ray emission is not associated with Liller 1. In fact 


Tam et al 


( 20111 ) localized the position of the intense y-ray source at ~ 7.5' from the cluster 


center. While this distance is smaller than the tidal radius quoted in H96 ( r t = 720" = 12'), 
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it is larger than the cluster tidal radius obtained in this work (r t ~ 300" = 5'). Even by 
taking into account the 7 -ray position error (~ 2' for the Fermi-LAT observations), the 
location of the 7 -ray source is only marginally compatible with the newly determined radial 
extension of Liller 1, thus opening the possibility that such a strong emission is not coming 
from this stellar system. 

This research is part of the project Cosmic-Lab (http://www.cosmic-lab.eu) funded by 
the European Research Council under contract ERC-2010-AdG-267675. D.G. gratefully 
acknowledges support from the Chilean BASAL Centro de Excelcncia en Astrofisica y 
Tecnologfas Afines (CATA) grant PFB-06/2007. F.M. gratefully acknowledges the support 
provided by Fondecyt for project 3140177. S.V. gratefully acknowledges the support 
provided by Fondecyt reg. 1130721. This work made use of data products from the Two 
Micron All Sky Survey, which is a joint project of the University of Massachusetts and the 
Infrared Processing and Analysis Center/California Institute of Technology, funded by the 
National Aeronautics and Space Administration and the National Science Foundation. 


Tabic 1: New parameters for Liller 1. 


Parameter 

Derived Value 

Center of gravity 

c«j 2 ooo = 17 h 33 m 24.56 s 


dj2ooo = —33°23'22.4" 

Reddening 

E(B-V) = 3.30 ±0.20 

Distance Modulus 

(m - M) o = 14.55 ± 0.15 

Distance 

cl = 8.1 ± 1.0 Kpc 

Core radius 

r c = 5.39" = 0.21 pc 

Effective radius 

r e = 22.4" = 0.88 pc 

Half-mass radius 

r h = 30.5" = 1.20 pc 

Tidal radius 

r t = 298" = H .74 pc 

Concentration 

r — 1 74 _015 

c — +0.15 

Total luminosity 

ibol = 6.9 +g x 10 5 Lq 

Total mass 

M tot = 2.3 +3 x 10 6 M 0 (Salpeter IMF) 
M tot = 1.5 ;°;2 x 10 6 M 0 (Kroupa IMF) 

Central mass density 

po ~ 7.2 x 1O 6 M 0 (Salpeter IMF) 

po ~ 4.8 -{% x 10 6 M 0 (Kroupa IMF) 

Central RT s -band SB 

Pk b { 0) = 10.29 ± 0.17 mag arcsec -2 

Central J-band SB 

p.j( 0) = 12.75 ± 0.17 mag arcsec -2 
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Fig. 1.— False-color image of Liller 1 obtained by combining GEMINI observations in the 
near-IR J and K s bands. North is up, east is on the right. The held of view is 85" x 85". 
The bluish stars are hot held objects. 
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SOFI (J-K.) 

Fig. 3.— Distributions of the (SOFI - GeMS) magnitudes as a function of the SOFI color 
(J — K s ). Stars at a distance larger than 30" from the cluster center and with 2.8 < (SOFI K s 
- GeMS K s ) < 3.6 and 2.6 < (SOFI J - GeMS J) < 3.6 respectively, are shown in figure. The 
solid lines instead represent the median values obtained by applying an iterative u-clipping 
procedure. Finally we obtain (SOFI K s - GeMS K s ) = 3.299 and (SOFI J - GeMS J) = 
3.146. 
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Fig. 4.— Near-IR differential reddening-corrected CMD of Liller 1 obtained from the GEM¬ 
INI observations discussed in the paper. The main evolutionary sequences of the cluster are 
well visible down to the MS-TO point. On the blue side of the CMD, for (J — K s ) < 1.5, the 
blue plume defined by the Galactic held MS is also distinguishable. The photometric errors 
for each bin of K s magnitudes are shown on the right side of the panel. 
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Fig. 5.— Near-IR CMDs of Lillcr 1 as obtained from the VVV data. In the left panel only 
stars external to the GEMINI FoV and located at distances smaller than r t from C grav are 
shown. In the right panel only stars located at r t < r < 500". 




Fig. 6.— Luminosity function of the red clump and the brightest portion of the RGB in 
Lillcr 1 (upper panel) and in the reference cluster NGC 6553 shifted by 5K = 1.9 (lower 
panel). The two dashed vertical lines mark the location of the red clump and the RGB 
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Fig. 7.— Star density profile of Liller 1, obtained from resolved star counts in the combined 
data set. Empty circles represent the observed profile, while solid ones are obtained after 
subtraction of the Galactic held density (marked with the dotted line). The best-fit single¬ 
mass King model is shown as a solid line and the corresponding structural parameters are 
labeled in the figure. The lower panel shows the residuals between the observations and the 
fitted profile at each radial coordinate. 
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Fig. 8.— SB profiles of Liller 1 in the K s and J bands (top and bottom panels, respectively), 
obtained from the 2MASS data set. Empty circles represent the observed profile, while solid 
ones are obtained after background subtraction (marked with a dotted line). In both panels, 
the solid line corresponds to the best-fit King model shown in Figure [3 
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